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An effective Hamiltonian describing the surface states of a toroidal topological insulator is ob¬ 
tained, and it is shown to support both bound-states and charged zero-modes. Actually, the spin 
connection induced by the toroidal curvature can be viewed as an position-dependent effective vec¬ 
tor potential, which ultimately yields the zero-modes whose wave-functions harmonically oscillate 
around the toroidal surface. In addition, two distinct Berry phases are predicted to take place by 
the virtue of the toroidal topology. 
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I. INTRODUCTION AND MOTIVATION 

Topological quantum states of matter do not break 
any local symmetry, so that they cannot be described 
by Landau theory of phase transitions. Instead, they 
demand topological quantum numbers associated with 
the bulk wavefunctions. Examples include the integer 
quantum Hall effect, topological insulators (H, topo¬ 
logical superconductors and topological superfluids j2|. 
On the other hand, geometrical features of the real and 
momentum spaces are very important to describe the 
physical properties of several condensed matter systems, 
like the nematic order in liquid crystals Q, graphene 


where R is the sphere radius EMI- Furthermore, the 
effect of the spin connection in surface carriers can be 
interpreted as fictious magnetic monopole and there 
are two types of Berry phase for the surface electrons 
0 - It is noteworthy to mention that the case of a 
spherical TI has been considered by fl3j ], while a more 
general situation, describing any curved topological 
insulator was presented in Ref. [16j. Here, we study 
the toroidal geometry, whose variable curvature leads 
to an induced spin connection which may be faced as 
a position-dependent effective vector potential. Fur¬ 
thermore, such a feature yields the zero-modes and two 
distinct Berry phases to surface carriers. In addition, a 
conical TI presents a electric polarization associated to 
its aperture angle 0 , so that narrower or wider cones 
appear geometrically polarized in opposite directions. 
Other studies include cylinder 0 and the Mobius strip 
[20j . the simplest non-orientable Riemann manifold, 
where the usual topological invariant is not well-defined. 
Moreover, the classical limit has been used to study the 
scattering of electrons in curved TI’s with non-trivial 
metrics (17j |. 

Despite the increasing interest in the study of curved 
TI’s, our knowledge is still very superficial with several 
questions demanding answers. For instance, one knows 
that the Dirac cones are gapless along all orientations 
of a flat TI surface [in]]. Whether this gapless feature 
is kept whenever an arbitrary curved TI is considered 
comes to be an open question, even though the sphere 
has provided a partial and affirmative answer, as pointed 
out above. At this point, the toroidal geometry appears 
to be a very suitable manifold to be studied, once it 
encompasses a variable curvature interpolating between 
negative (pseudosphere or hyperbolic plane) and posi¬ 
tive (sphere) curvatures. In turn, such a geometry has 
recently received a considerable attention in the litera- 

_ ture. As examples, we may quote toroidal tight traps 

for Bose-Einstein condensates (22j, carbon nanotubes 
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Due to their potentialities for application in spin- 
tronics and quantum computation, there has been an 
increasing interest in the study of topological insulators 
(TI’s) in the last years. TI’s exhibit an insulating bulk 
while their borders support conducting states protected 
by time-reversal symmetry (TRS) [lj. A TI is character¬ 
ized by a topological number associated with the Bloch 
wavefunctions describing the bulk band structure which 
can assume two values, 0 or 1, accounting for an ordinary 
or a topological insulator respectively. In fact, a TI has 
an odd number of Dirac cones connecting the electronic 
energy bands on the TI surface and may be described by 
a two-dimensional Dirac equation ll,|12j. These cones 
lie over all orientations of flat surfaces, making them 
gapless. 

Actually, curvature has been shown to affect the 
physical properties of topological insulators. For in¬ 
stance, on a spherical TI, it has been reported that the 
surface carriers remain gapless while the bound-states 
experience a gap proportional to the curvature, 1 /R, 




2 


can support a vortex as the magnetization ground state 
for smaller radius than cylindrical nanorings [2(1 [27j ■ 

Our article is outline as follows: in Section II we 
present the model and obtain the surface effective 
Hamiltonian for a toroidal TI. Section III is devoted to 
study how the toroidal geometry and topology leads to 
two distinct Berry phases and the fictitious magnetic 
monopole induced by the curvature. Next, in Section 
IV, we discuss upon the physical spectrum, namely, we 
show that there is charged zero-modes, a fact brought 
about by the variable curvature of the toroidal surface. 
Finally, the conclusions and prospects are presented. 
In the appendix we present an alternative derivation 
of the surface effective Hamiltonian for a toroidal TI 
(supplementary material of section II). 


II. SURFACE EFFECTIVE HAMILTONIAN FOR 
TOROIDAL TOPOLOGICAL INSULATOR 


We consider BHZ model describing TI’s such as Bi 2 Se 3 , 
Bi 2 Te 3 and Sb 2 Te 3 |l!,[29j, with a minimal set of physical 
parameters which encompasses the distinction between 
an ordinary and a topological insulator. This is done 
by requiring particle-hole symmetry and isotropic mass 
term, as below: 


U{k) 


M(k) Ak z 0 Ak- 
Ak z —M(k ) Ak _ 0 

0 Ak+ M(k ) — Ak z 
_ Ak + 0 —Ak z —M{k). 


(1) 


where k± = k x ± ik y and M(k) = M — B(k% + k y + fcf); 
M/B < 0 concerns for an ordinary while M/B > 0 ac¬ 
counts for a topological insulator. It is noteworthy that 
even in this simplest case, the surface states appear when¬ 
ever there is a topological phase transition, say, whenever 
M/B changes its sign at the boundary, for example, pos¬ 
itive in the bulk and negative in the vacuum. These 
parameters are computed by fitting the energy spectrum 
of the effective Hamiltonian to that coming from ab ini¬ 
tio calculations [28j. In Table [0 we quote typical values 
of such parameters for a number of common topological 
insulators. 



FIG. 1. [Color online] Above, a ring torus embeded in 3D 
space is depicted; Below, a cut showing the peripolar coordi¬ 
nates, (9, (f)) € [0, 27t), and the internal and external torus 
radii, r and R. 


(a single central hole), and whenever embedded in three 
dimensional space, it shapes like a donut (see Fig. [I]). 
The standard tori are classified in three distinct types 
concerning the relations between their internal, r, and 
external, R , radii. For R > r, we have a ring torus 
(donut shape, as shown in Fig. QJ, if R = r a horn torus 
is obtained, while for R < r a spindle torus takes place. 
Any ordinary torus may be parametrized, for instance, 
by peripolar coordinates ( 9 , <j>): 


(R - \jx 1 + y 2 ) 2 + z 2 =r 2 , (2) 

which are related to the Cartesian ones by: 



A (eV A) 

B (eV A 2 ) 

M (eV A) 

Bi 2 Se 3 

2.26 ~ 3.33 

6.86 ~ 44.5 

0.28 

Bi 2 Te3 

0.30 ~ 2.87 

2.79 ~ 57.38 

0.30 

Sb2Te3 

0.84 ~ 3.40 

19.64 - 48.51 

0.22 


TABLE I. Typical values of physical parameters A and B for 
some TI crystals in the topological insulator phase [I]. 

Our interest is to study the above model on the toroidal 
geometry, which is a smooth surface with varying curva¬ 
ture. The simplest torus is a surface having genus one 


{ x = (R + r sin 9) cos (f> , 

y = (R + r sin 9) sin cj) , (3) 

z = r cos 9 , 

where R is the rotating (external) and r is the axial (in¬ 
ternal) radii, respectively (see Fig[T|). Thus the momen¬ 
tum operator in Eq. m is written in toroidal coordinates 
as: p = p r r + pgO + where p r = —id r , pg = —idg/r 
and p</, = —id'frKR + r sin 6). 

In what follows, we shall derive a surface Hamilto¬ 
nian in the k ■ p approximation. First, note that T~i(k) 
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may be split into two parts: 7i(k) = T-L\\ + 'H±, where 

= R(k) | P9= o,p^,=o- Let us begin, by solving the ra¬ 
dial eigenvalue problem: 


The effective surface Hamiltonian acts onto the spinor 
| a >, which in the k ■ p approximation, may be written 
as |a) = a-i-|+} + a_|—}, so that 


H.lIV’) = E^) . 


Assuming that | ip) may be written as: 


\iP) = \iP(r,9,cp)) = e x ^\u(9,cp)), 


( 4 ) 

( 5 ) 


" <+ 73|||a)" 

= 'H s \a> with H s = A 

0 i>+ 



£>_ 0 


( 11 ) 


whi 


where A accounts for the penetration of the surface wave 
function into the bulk. Demanding that such a solution 
satisfies the boundary condition, we must have: 

\iP(r = R)) = 0, (6) 

which states that all four components of the wave func¬ 
tion vanish on the surface of the torus (at r = R). 
This boundary condition is analogous to that used in the 
spherical caseJUJ ensuring that there are no charge car¬ 
riers outside the TI surface. Boundary condition above 
demands: 


A± 


-a/3 A ± sjA 2 + 4MB 
2 B 


( 7 ) 


where (a,/3) = ±1 specify the spin polarization of the 
radial eigenvector. Since ip describes localized surface 
states, A± must be both positive, which demands a/3 < 0 
and MB > 0, since A/B > 0 (see Table HI). In addition, 
note that the requirement M/B > 0, coming from 
boundary conditions, must be fulfilled in the topological 
state [l|. 


Solution of Eq. © may be written as (N is a normal¬ 
ization constant): 

ip = N[e x + {r ~ R) - e A -( r -- R )]u[A + ] = p(r)u[A+]. (8) 


Now, by defining the basis eigenstates of H±, as \ip >= 
|±), where 


|+) = p(r)|f+}_ |-) = p(r)|r-)+, (9) 


we obtain: 


|f+}_ = | 


1 

—i 

1 

—i 


e cos | 


e'sin § 


IM+ = 5 


e“^/ 2 sin | 


e ^/ 2 cos | 


( 10 ) 


- dg 1 cos 9 i 

7A± —+ T „ „ . Q + D . 

r 2 R + r sm 0 R + r sin 0 

which is the Dirac operator for a free massless fermion on 
a toroidal surface. It is noteworthy that whenever R —> 0 
we recover its spherical counterpart [13]. Instead of using 
the double-valued basis we opted by a single-valued 

one, then we should replace —idj by —idj + 1/2 (j = 6,<p). 

III. BERRY PHASES IN A TOROIDAL TI 

The two distinct Berry phases appearing on the surface 
of a toroidal TI are associated to the angles 6 and <p. 
Indeed, this comes about from the fact that a torus may 
be generated by the topological product of two circles, 
say, Sg # j30|. To see this in more details, let us note 

that any solution of 

H s ip{0, (p) = Eip{6 , <p ), (13) 

must satisfy periodic boundary conditions: 

ip(9, <p) = ip(9 + 27t, cp) = ip(9, <p + 27r). (14) 

Here, ip — I 0 ) = a +l+) + a_|—}, and using Eqs. (JSJ- 
HQD, we obtain: 

a±(9, (p + 27r) = a±(9 + 27r, <p) = —a±(9 , <p). (15) 

These anti-periodic boundary conditions are equivalent 
to Berry phases, 7 = 7 r, as follows: 

a{6 + 27t, <p) = e lle a{6 , <p) = e l7r a(9 1 <p) 
a( 0 , (p + 27t) = e* 7 ^a( 0 , cp) = e i 7 r a(d, <p). (16) 

Whenever a charge carrier, whose spin is locked to the 
surface, completes a closed path along the azimuthal or 
the polar angle, any of such loops cannot be shrunk to a 
point by means of smooth deformations. Physically, its 
spin is twisted by 2n in such a way that ip accumulates a 
7 r-phase. Indeed, TI’s have a nontrivial intrinsic momen¬ 
tum space Berry phase acquired by an electron when it 
revolves around a Fermi loop for its spin rotates with the 
momentum k around Fermi surface. This is important to 
describe the electron behavior. For instance, whenever it 
couples to a magnetic field or it is subject to disorder ef¬ 
fects, the 7r-phase ensures antilocalization, that is, even a 
strong disorder cannot localize surface electronic states, 
provided that the bulk energy spectrum remains gapped. 
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Alternatively, the Berry phases may be related to the 
spin connection, Tg = which is equivalent to 

a vector potential, A = Tgcf>, generated by an effective 
magnetic charge distribution induced by the toroidal cur¬ 
vature, K = sin (9)/r(R + rsin(0)). Such a distribution 
is expected to be a ring of magnetic monopoles, with 
radius R, placed on the xy-plane. If we take R —> 0, 
it shrinks to a single Dirac magnetic monopole, located 
at the sphere center, as obtained in Ref. 13]. In addi¬ 
tion, the curvature appears as a geometrical mechanism 
to induce artificial gauge fields in TI as much as mechan¬ 
ical tension and topological defects do in the graphene 
framework |3j]. 


IV. PHYSICAL SPECTRUM: BOUND-STATES 
AND ZERO-MODES 

The non-trivial curvature of a manifold brings about 
spin connection on the Dirac equation. Namely, this af¬ 
fects the electronic properties of surface carriers on a 
toroidal geometry: energy spectrum, spin polarization, 
local density of states, and so forth. Let us consider the 
eingenvalue problem m with the following ansatz 


a(0,$=e im * ( am + f Q \ ) , (17) 

\ a m -(0) ) 

where m accounts for the azimuthal component of the or¬ 
bital angular momentum, L§, taking half-integer values, 
accordingly to boundary conditions m , say, to = ±N /2 
with N an odd integer. Thus, in terms of a m ±(6) equa¬ 
tion above becomes 


Id 1 cos# 

r dd 2 R + r sin 6 


TO 

R + r sin#. 


— 


E 

~7 a m + ■ 
(18) 


Id 1 cos# to 1 E 

r dd 2 R + r sin 9 R + r sin 9. a ™ + A 

(19) 

These differential equations may be decoupled at second 
order, as follow: 


1 d . , . d sin# 

_ (p + sin 9) d9 ^ ^ rf# 2(p + sin#) 


1 


(p + sin#) 2 


TO + — COS d ) + £ 


= 0, (20) 


where a = ±, e = rE/A, and p = R/r. To our best 
knowledge, this differential equation presents no known 
analytical solution for e ^ 0 (zero-modes, e = 0, will 
be discussed below). Efforts have been done to solve it 
numerically; however, such findings have not been shown 
useful for a better comprehension of the bound-state 


spectrum. On the other hand, a qualitative analyses 
indicates that the energy gap depends on the internal r 
and external R radii of the torus and it is characterized 
by two independent quantum numbers, to and n , whose 
discrete energy levels E(n,m ) must be symmetrical 
around E = 0 (by virtue of the electron-hole symmetric 
theory). In addition, let us note that as p = R/r —> 0 the 
spherical surface eigenstates along with their spectrum 
are recovered [l3|, as expected . 

Zero-modes: Although surface states on the torus 
could not be solved analytically for an arbitrary energy, 
zero-modes solutions may be readily worked out. Actu¬ 
ally, by taking E = 0 in eqs. m and m, we clearly 
realize that they are solved by: 

a 0 (9,<f,)=Ne im *( am + ( °\ ) , (21) 

V “m-W J 


Qm + iP ) — 


r r?_ 


(xu-S+A 


sin 9 


, ( 22 ) 


a m - (9) = \/ ^ + sin 9 e 


^ V^- 1 ^ . (23) 


each of them bearing an integer electronic charge: 

Q/e~ = J d 2 rj° = J Gprajao = 1. (24) 


The zero-modes states appear with wave-functions that 
harmonically oscillate around the toroidal surface in the 
azimuthal and polar directions ((f) and 9 coordinates). To 
explain the physical origin of these zero-modes we recall 
that the toroidal curvature is equivalent to an effective 
vector potential, A = f ^_)]° s s f n g <t>, as pointed out in the 
previous section. Now, let us consider a given zero-mode 
displacing a loop along <f >, in such a way that it expe¬ 
riences an effective magnetic flux given by (recall that 
dl = drr + rd9 9 + (R + r sin 9)d<f> <j >): 


$ = A ■ dl = 7T cos 9. 

Therefore, a closed path along <f> with constant curva¬ 
ture, one realizes a zero-mode state, labeled by a quan¬ 
tum number to. The curvature then plays a role analo¬ 
gous to an external field axially applied to a TI cylinder, 
where it is induced a pair of zero-modes depending on 
the magnetic flux [32:]. Furthermore, once the torus cur¬ 
vature smoothly varies as 9 goes from 0 to 27r, it offers 
a unique framework where zero-modes appear due to the 
geometry of the surface. 

For the sake of completeness, it is worthy to men¬ 
tion that the emergence of zero-energy modes can be 




























5 


established, a priori, by mathematical index theorems, 
which relates the occurrence of such modes to the ge¬ 
ometry and topology of the space on which the Dirac 
equation is stated jUj, HU- For instance, a topologically 
non-trivial background or a defect may bring about a 
position-dependent mass term into the Dirac equation 
enabling the spectrum to support isolated zero-modes 
[35l |36j. Charged zero-modes can be also found in a 
topological insulator when its surface is coated with a 
ferromagnetic film supporting an in-plane magnetic vor¬ 
tex. In this case, they have the additional property of 
being fully polarized modes located near the vortex cen¬ 
ter |33]. 


V. CONCLUSIONS 

We have shown that the surface of a toroidal topo¬ 
logical insulator supports zero-energy modes, while the 
bound-states have an energy quantization depending on 
two independent quantum numbers, each one associated 
with a distinct circle composing the torus by means of a 
2-circle topological product, T 2 = S 1 ffS 1 . Actually, the 
torus variable curvature leads to an additional contribu¬ 
tion to the surface states eigenvalue problem, preventing 
us from obtaining analytical solutions with non-vanishing 
energy. 

On the other hand, we have realized that it is precisely 
this variable curvature which enables the zero-energy 
spectrum to appear, as the case of a TI cylinder subject 
to an applied magnetic field [32;]. In our framework, the 
curvature mimics the effect of a vector potential, yield¬ 
ing non-trivial effective magnetic flux to the surface elec¬ 
trons whenever they displace along a closed loop along 
the torus major radius, say, angle 4>. Therefore, for each 
closed path with a fixed curvature, the modes evolve as if 
they were experiencing a given flux. Therefore, our find¬ 
ings suggest geometrically non-trivial manifolds as useful 
frameworks for the study of how the curvature (and per¬ 
haps other ingredients, like torsion) affects the physical 
properties of surface states in topological insulators. 

In practice, the best fabricated samples are still imper¬ 
fect in geometry, pureness, and so forth. Even for a re¬ 
alistic toroidal TI, where (small) imperfections and non¬ 
magnetic impurity concentration take place, our main 
results are not expected to be jeopardized. For instance, 
the zero-modes surface states must be kept, once it relies 
on more general considerations, say, the torus topology 
and its smoothly variable curvature. 

As a final remark, it should be stressed that it is pos¬ 
sible to construct a tight-binding model, taking into ac¬ 
count nearest-neighbor interaction, as it has been done 
for the spherical case. By solving such a Hamiltonian, its 
physical spectrum may be employed even for very small 
systems, say, nanoscaled torii, presumably with the sim¬ 
ilar properties as found here. 
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Appendix A: Dirac equation on a torus 

In order to confirm the results presented in Section II, 
we shall obtain the Dirac equation for a toroidal TI fol¬ 
lowing the general method from Ref. (l6| . There, it is re¬ 
alized that in a curved surface the Dirac equation is modi¬ 
fied by the spin connection. On the other hand, a genus-1 
torus, Fig. [lj is a compact surface whose Gaussian curva¬ 
ture is given by K = , smoothly varying from 

— l/r(R— r) to +l/r(R + r) along its polar (0) angle, so 
interpolating between the pseudospherical and spherical 
curvatures whenever R > r, and vanishing when 9 = 0 
or 9 = ir. Despite the Dirac equation on a torus can be 
written using the covariant formalism of quantum field 
theory in curved space time, here we adopt the general 
approach of Ref. 16] which is equivalent to the method 
presented in the previous sections. In Ref. [jj| the effec¬ 
tive Hamiltonian for the two-component spinor a (0, <j>) 
is given by: 


Us 


0 V+ 
£>_ 0 


(Al) 


where 

2 

V+ = Y, [(°‘- A ~ b ^ B ) ( d * + 2 dl ln ^)) 

2=1 

+ 2 l d i( a i A ~biB) , 


(A2) 


£>_ = 


^ 1 
X] [ - - W ( d * + 2 di ln ^>) 


2=1 


- -[di(a,iA- biB)* 


(A3) 


The parameters and bi are determined by the torus 
geometry and are given by: 


1 

a\ = — , a 2 
r 


i 

R + r sin 9 ' 


bl ~ ' 


b 2 — — - 


i sin 0 


Acknowledgments 


(.R + r sin0) 2 
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Assuming that the penetration depth A for the surface 
states is much shorter than r, then (r) can be approxi¬ 
mated by the surface radius and (VG) = r(R + rsin0) 
(see Ref. [n| for more details). Therefore, Dirac opera¬ 
tor V± reads: 



dg 1 cos 9 
^ r R + r sin 9 


i(A — B sin#) 
R + r sin 9 


(A4) 


where a^A — biB may be faced as an effective velocity in 
x l = (9, 4 >) direction; B/r is a (renormalization) correc¬ 
tion to the effective velocity of the surface carriers, while 


B sin 9 accounts for the renormalization due to curvature. 
However, these terms can be ignored at first approxima¬ 
tion, in such a way that Dirac operator (1A4I) comes to 
be identical to the surface Hamiltonian m, derived in 
Section II. 

It should be emphasized that a factor = r(R + 
rsin#), (Q = det {gtj}, gij is the torus metric) must be 
introduced in the definition of the matrix elements of 
the Hamiltonian GD whenever evaluated between sur¬ 
face states, because in a 2D manifold the natural integral 
measure is dS = \Z r Qdx 1 dx 2 = r(R + r sin 9)d9d(j). The 
correction term | R+° s s e m g hr the above linear differential 
operator corresponds to a spin connection in the Dirac 
theory on curved surfaces [l6j, describing the effect (cou¬ 
pling) of the curvature on electrons motion. 
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